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The Duru- Kleinert (DK) method of solving unknown path 
integrals of quantum mechanical systems by relating them 
to known ones does not apply to Markov processes since the 
DK-transform of a Fokker-Planck equation is in general not a 
Fokker-Planck equation. In this note, we present a significant 
modification of the method, based again on a combination of 
path-dependent time and coordinate transformations, to ob- 
tain such relations after all. As an application we express un- 
known Green functions for a one-parameter family of Markov 
processes in terms of the known one for the Schenzle-Brand 
process. 

1. The stochastic theory of Markov processes ex- 
plains many phenomena where fluctuations play a signif- 
icant role. Prominent examples are provided by the emer- 
gence of self-organization ||, the occurrence of quan- 
tum dissipation 0] and the appearance of stochastic res- 
onance [^|. In accordance with their wide range of appli- 
cability there exist various powerful solution methods for 
Markov processes based either on the global characteriza- 
tion of the probability evolution by the Onsager-Machlup 
path integral or on its local equivalent, the Fokker-Planck 
equation. Examples for the former are the small noise ex- 
pansion and the adiabatic elimination procedure of fast 
random variables, for the latter the eigenfunction expan- 
sion and the continued-fraction method applied to peri- 
odically driven systems 

In this note we demonstrate the use of path-dependent 
time transformations, which have become a powerful tool 
for solving quantum mechanical problems since Duru and 
Kleinert 's original work on the path integral of the hy- 
drogen atom wjfl. In one dimension with coordinate q, 
the crucial time transformation has the form 



as 



(1) 



where f(q) is some positive but otherwise arbitrary func- 
tion. Such a transformation does not change the standard 
formulations of quantum mechanics, a form invariance 
which has recently been emphasized in [p|JTo| . Due to 
this form invariance different quantum mechanical sys- 
tems can be related to each other. 

Some years ago, the DK-method was also applied to 
the stochastic theory of Markov processes 0,^2) , thereby 
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relating the Fokker-Planck equation of Markov processes 
to other stochastic differential equations. The latter had, 
however, an important disadvantage: they were no longer 
Fokker-Planck equations so that the DK-transformations 
did not link different Markov processes. This defect will 
be eliminated in the sequel by a significant modification 
of the DK-method. 

2. Consider a one-dimensional Markov process of a single 
random variable x, whose conditional probability density 
P{x 1 XQ]t) possesses the initial condition 

P{x,x o ;0) = S(x-x ) (2) 

and obeys the Fokker-Planck equation for t > 

d t P(x,x ;t) = H{x)P(x,x ;t). (3) 

Here H (x) denotes the infinitesimal time evolution oper- 
ator 



(4) 



containing the drift coefficient K(x) and the diffusion 
coefficient D(x). Then the Laplace transform 



G(x, xq; E) = / <il( 




-El 



P(x,x ;t) 



(5) 



represents a fixed-energy Green function solving the 
time-independent equation 



E 



(i) 



G^(x,x ;E^) + S(x-x ) = 0. (6) 



We have added a superscript (i) to emphasize that this 
equation describes the initial stochastic system from 
which we depart. 

It is well-known ||] that the Fokker-Planck equation 
remains form invariant under arbitrary invertible coordi- 
nate transformations 

X = x(q) , (7) 

where the Green function transforms as 

G^(q,qo;E^)=±x'(q)G^(x(q),x(q );E^). (8) 

The different signs take into account whether ^ is 
monotonously in- or decreasing. In fact, we conclude 
from (|) and'(§) that (q, q ; E®) satisfies the Fokker- 
Planck equation 



1 



ffW(g) - E®] GM(q,qo;E®) + 5(q - q a ) = 0, (9) 



where the new infinitesimal time evolution operator 



d 



tfW(g). = -— K {1 \q) 



1 d 2 



2 8q 2 

contains the drift and diffusion coefficients 



D^(q).] (10) 



^ (1) W = ^ (i) (^))-^|y^(^)), (ID 



D^{q) = -^D^{x{q)). 



(0/ 



(12) 



These coordinate transformations are of standard use in 
finding unknown solutions from known ones [ |l3||l4| ] . 

Let us now supplement these transformations by the 
path-dependent time transformation (|l|). First we pro- 
ceed in analogy with [|[ Chap. 12], and change the Green 
function according to 



G^(q,q ;E^) = f(q) 



F(q ;EW) 



xG^(q,q ;E^(E^)), 



(13) 



where F(q;E®) and E^(E^ ) are as yet unknown trial 
functions. Applying (0) with (|To|), we find the equation 
for the final Green function 



H^(q)-E ( - i \E^)+X(q- : E^)\ G«(g,g ;i; (f) (£ (i) )) 
+ <J(g-g ) = 0, (14) 

with the infinitesimal time evolution operator 



H {l Hq). = -^[K^(q), 



2cV 



£> (f) (g)« 



(15) 



containing the transformed drift and diffusion coefficients 



K (i \q)=f{q) 
D^(q) = f(q)D^(q) 
and the additional term 



F(q;E®) 



(16) 
(17) 



X(g;£«)=/(g) 



2 W F(q-EM) 



where 



f(q) 



(18) 



(19) 



The equation (|1J) has the above-mentioned defect of 
not being a Fokker-Planck equation, due to the presence 



of the additional term X(q;E®). This term can, how- 
ever, be removed by choosing any functions i?W(2?W) 
and F(q; E^ 1 ') which solve the differential equation 



(20) 



Note that although this equation is of the same complex- 
ity as the initial Fokker-Planck equation, only a partic- 
ular solution is required, so that labor will definitely be 
saved by our method. 

An alternative procedure which avoids solving the dif- 
ferential equation ( p0| ) is by leaving the time transfor- 
mation function /(g) open, choosing some trial func- 
tion F(q;E®), and calculating £(q;E®) from Jp|) with 
(pp|). If this happens to be factorizable as in (|19|), the g- 
dependent prefactor may be chosen as the transformation 
function /(g) and the E® -dependent one as the energy 
function EW(E®). 

In either procedure, the function F(q;E®) is subject 
to an important restriction. In the limit E® — > it has 
to satisfy 



lim F(q;E (i) ) = 1 



(21) 



identically in q, so that the energy function E^(E^) 
obeys 



lim eW(E®) 

E«->0 



0. 



(22) 



Only under this condition do initial and final Green func- 
tions possess proper stationary limits 



p£(x)= lim E®G®(x,x ;E®). 



EV>->0 



p$(q)= lim £( f )G( f )( Mo ;£ (f) ) 



(23) 
(24) 



/From (|), (H) and ©-© we read off a relation be- 
tween them 



J$C«0=± 



dE^jE® 



4^p£>fa(«))> (25) 



B(D=0 a: '(9( a; )) 

which guarantees the normalization of the probability: 
p«(x)dx= / p$(q)dq = l. (26) 



An interesting feature of the present method is that it 
permits us in the stationary limit to relate the proba- 
bility distributions of two arbitrary Markov processes to 
each other. Given initial and final drift and diffusion 
coefficients K(x) and D(x), we satisfy ( |l2| ) and ( |l7| ) by 
choosing the time transformation function as 



^ (f) (g) 
Dd)(x(q)) 



(27) 



f{q)=x u (q) 



Using this together with (y), fll2J), Ql6|), (|17j) and (|2lj), 
we obtain the desired coordinate transformation from the 
differential equation 



2 



x (q) = C exp 



x(q) 



1 „2K^(x) 

ClX — TTT ; r 



_2JTW(g) 



D(f)(g) 



(28) 



where C is an integration constant. The final stationary 
solution 



jy( f ) 

D(f%) 



exp 



D(f)(g) 



(29) 



is then related to the initial one 



■ exp 



, 2K®(x) 

dX ... ; ; 



(30) 



by (|25|), (|27D and (|28|), if the normalization constants 
satisfy 



C 



N( { ) dE® 



(31) 



BW=Q 



3. In order to demonstrate the applicability of the new 
transformation method we consider a Markov process 
with a multiplicative noise for a random variable x £ 
(0, oo ) where the drift and the diffusion coefficient de- 
pend on an arbitrary parameter a > as follows: 

K^(x)=a^x-b®x 2a+1 , D^(x)^Q^x 2a + 2 . (32) 

Performing a transformation of the random variable (Q) 
and the time (|f|) with 

x(q) = / , /(g) - <Z 7 , (33) 

the associated differential equation (^0|) with ([L8]) is 
solved by a function 

F( g ;£«)=^ B(i) ), (34) 

if the parameters (3 and 7 are related according to 

2o^ + 7 = 0. (35) 

The function 6(E®) and the energy relation E^> = 
E^(E^ ) are determined by 



f(23«) = -^ 



6(0 QW 



a« 2a« 



(36) 
. (37) 



and (37) satisfy the correct limits 



Note that 
@ and (|| 

The transformed drift and diffusion coefficients then 
follow from @, @, @ and (jfjj) as 

tf<%) = a < f >g - ^V 2 "^ 1 , Z?( f )(g)=Q( f ) (Z 2 , (38) 



where 



6< f > 



6(0 

(i) 



/3 



0-1 



2/3 



a 

/3 2 



(39) 
(40) 
(41) 



These relations supply us with solutions of the Fokker- 
Planck equation for the one-parameter family of Markov 
processes (^) if we specialize 



a 



(42) 



Then the final Markov process (|38|) for a random variable 
q £ (0, 00) coincides with the well-understood Schenzle- 
Brand process [0], which is a standard model in non- 
linear optics and chemical reaction dynamics. It can be 
derived as an approximation to a number of different pro- 
cesses by adiabatically eliminating fast random variables 
in the limit of large external fluctuations. For instance, 
the Schenzle-Brand process describes the electrical field 
near a laser threshold, the multiplicative noise being due 
to inversion fluctuations. 

With standard methods ||, the Fokker-Planck equa- 
tion of the Schenzle-Brand process with (|3^) and (Jl2|) can 
be transformed to the Schrodinger equation of the Morse 
oscillator. As the quantum mechanical Green function of 
this system has been explicitly calculated from path in- 
tegrals [ p^|Jl6| , the Green function of the Schenzle-Brand 
process is known: 



G^(q, qo ;E^) 



v ( h^-L 4-1 & 

\ 2 + 4 2Q 



2Q(f) 



aV) 1 



exp 



b (f) T (l + fc(f)) 
b {{ Hql-q 2 ) 



2Q( f ) 



fb^q 2 



4 + 2Q<f)' 2 



©(g - go) + (g <-> go) 



where denotes the abbreviation 




2E(f) 



(43) 



(44) 



With the help of the transformation formulas (||), (|l|), 
and taking into account (]33|)~(f4"2"|), we find from this the 
unknown Green function of the initial Markov processes 



G^(x lXo ;E®) = 



-ak^ —a — 1 



XX, 



ak^ 1 ' +c 



exp 



a«r (1 + jfeW) 

aW / 1 1 



2aQ« 



3 



x{e(x Q -x)w k 



with 



n(0 



i) ^2a 



(45) 



feP) 



1 

2a' 



aQP) 

Furthermore we obtain from (M) 



(46) 

the station- 



ary solution of the Schenzle-Brand process ( |38[ ) , fl42J ) 



»(*> 

Q< f ) ' 



r f a (f) 

1 I Q(f) 



exp 



6 (f) 



(47) 



which is mapped via ( p5~| ) 
process (p2|), yielding 



to the initial Markov 



2a 



Q(i) 



-+1 



x.r 



r 

'q(') 



-2a-2 



exp 



,.P) 



(48) 



The analytic properties of the Green functions 
GP)(x,a;o;£P)) and G^(q, q ; E^) in the energies E® 
and 23P) determine the spectra of the infinitesimal time 
evolution operators H^(x) and H^ { \q) |8[ Chap. 9]. 
From (^3|)-( ^5|) we deduce that the initial multiplica- 
tive process ( |32|) has only a discrete spectrum, whereas 
the final one ([38[), ( (42| ) contains both a discrete and a 
continuous branch. Such differences between spectral 
types were encountered before in quantum mechanical 
DK-transformations: the hydrogen atom has discrete and 
continuous states whereas the DK-equivalent oscillator 
has only discrete states J?],!)]. The spectra are usually 
related by a Sommerfeld- Watson transformation of the 
Green functions §, Chap. 14], @. 

The discrete levels closest to zero are in a one-to-one 
correspondence HJll|. For stochastic systems, these lev- 
els rule the approach of the conditional probability den- 
sity to its stationary limit. In our example the associated 
poles of the initial Green function ( [45| ) 

£P } = -2aa^n , n = 0, 1, . . . (49) 

are mapped to the corresponding final ones of (|43|), fl44| ) 

=2Q( f V + (QP>-2aP))n (50) 

by 

£P) = (51) 
as long as n is bounded by 



n < 



qP) 
2QP) 



(52) 



The evaluation of (|Tj) requires the use of relations (|3 
and ©-©. 

4. In summary, we have shown that a combination of lo- 
cal time and coordinate transformations opens new possi- 
bilities of relating different Markov processes. By extend- 
ing the method to several random variables, we expect 
many useful applications. Furthermore we hope that this 
method might help to solve non-Markov processes |19|. 

Let us finally mention that the local time transforma- 
tion (Q) is nonholonomic in spacetime, i.e., it carries a 
flat spacetime into a spacetime with nonzero torsion and 
curvature |^,|l0|- If we allow for purely spatial nonholo- 
nomic changes of coordinates, we can reach also space- 
time geometries with spatial curvature and torsion. This 
will enable us to describe technically relevant diffusion 
processes in crystals with defects [p0|-p3|. 
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